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ABSTRACT 

Aims. We report the first detection of the J = 1 - 0 (102.6 GHz) rotational lines of CF^ (fluoromethylidynium ion) towards CygX-N63, 
a young and massive proto star of the Cygnus X region. 

Methods. This detection occurred as part of an unbiased spectral survey of this object in the 0.8-3 mm range, performed with 
the IRAM 30m telescope. The data were analyzed using a local thermodynamical equilibrium model (LTE model) and a population 
diagram in order to derive the column density. 

Results. The line velocity (-4 km s"^) and line width (1.6 km s"^) indicate an origin from the collapsing envelope of the protostar. 
We obtain a CF+ column density of 4xl0^^cm“^. The CF^ ion is thought to be a good tracer for and assuming a ratio of 10"^ for 
CF^/C^, we derive a total number of of 1.2x10^^ within the beam. There is no evidence of carbon ionization caused by an exterior 
source of UV photons suggesting that the protostar itself is the source of ionization. Ionization from the protostellar photosphere is 
not efficient enough. In contrast. X-ray ionization from the accretion shock(s) and UV ionization from outflow shocks could provide 
a large enough ionizing power to explain our CF^ detection. 

Conclusions. Surprisingly, CF+ has been detected towards a cold, massive protostar with no sign of an external photon dissociation 
region (PDR), which means that the only possibility is the existence of a significant inner source of C^. This is an important result 
that opens interesting perspectives to study the early development of ionized regions and to approach the issue of the evolution of the 
inner regions of collapsing envelopes of massive protostars. The existence of high energy radiations early in the evolution of massive 
protostars also has important implications for chemical evolution of dense collapsing gas and could trigger peculiar chemistry and 
early formation of a hot core. 
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1. Introduction 


The formation of massive stars is still not well understood. Im¬ 
portant questions are related to accretion rates required to form 
massive stars in the monolithic collapse scenario, and to the evo¬ 
lution of the inner regions of the collapsing cores where the 


disks form and stellar feedback o ccurs (e.g. |Yorke & Sonnhalter 
2002 Duarte-Cabral et al.||2013 ). For high accretion rates (10 
to 10 Mo/yr), the protostars have large radii (10 to 100 R©) be¬ 
fore the y contract to reach the main sequence while still accret¬ 
ing (e.g. |Hosokawa & Omukai|2()09] ). It is roughly when the pro¬ 
tostars reach the main sequence that they start to develop ultra¬ 
compact (UC HII) regions which may then regulate and possi¬ 
bly stop further accretion, limiting the final mass of the stars. 
Massive protostars also form hot cores with large fractions of 
collapsing gas radiatively warmed over 100 K. They release into 
the gas phase molecular species previously frozen on grains and 
drive a warm gas chemistry. It is not clear whether this hot core 
phase occurs before or at the same time as the development of 
an UC HII region. By opening large inner cavities, outflows and 
UC HII regions may help to heat up a large fraction of collapsing 
envelopes, leading to well-developed hot cores. 


CygX-N63 has been determined by |Motte et 11] ( [2007| ) to 
be one of the brightest and most compact 1.2mm cores in the 
Cygnus X region ( [Schneider et al.|2006] ) at a distance of 1.4 kpc 
( jRygl et al.|2012| ). It is very likely a massive p rotostar in its earli¬ 
est phase of formation ( jBontemps et al.|20T0| ). With an envelope 
mass of 44 M© within 2500 AU, and a luminosity of 340 L©, it is 
the most massive and youngest protostar detected in Cygnus X, 
and is recognized as a Class 0 massive protostar driving a power¬ 
ful CO outflow ( jPuarte-Cabral et al.|2013| ). It does not yet excite 


an UC HII region and does not contain a developed hot core and 
could thus be a rare example of a massive protostar in its pre- 
UC HII region and pre-hot core phase. This exceptional object is 
one of the best known examples of a dense core in monolithic 
collapse with no sign of sub-fragmentation from 0.1 pc down to 
500 AU, although it could still form a close binary. 

Molecular lines are important probes of the dense gas from 
which stars form and hence provide important insight for under¬ 
standing the earliest phases of star formation. One of the major 
cooling lines of the warm interstellar medium is the far-infrared 
flne stru cture line of ionized carbon (C^) at 158 pm. It was rec¬ 
ognized ( jNeufeld et al.|2005] ) that in the molecular gas phase of 
the interstellar medium, HE is the main fluor reservoir and that 
CF^ could be a good indirect tracer of C^, which is not reach- 
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able from ground-based facilities. The first (and o nly) detections 
of mm-transitions of CF^ towards the Orion Bar dNeufeld et al 


|2006| ) and the Horsehead nebula (Guzman et al. 2012) confirmed 
these predictions and indicate that CF^ can trace in photon 
dominated regions (PDRs). 

Here we present the first detection of CF^ towards a protostar 
obtained within a 0.8 to 3 mm unbiased spectral survey employ¬ 
ing the IRAM 30m towards CygX-N63. In Sect. 2 the observa¬ 
tions are detailed while results are provided in Sect. 3. Sect. 4 
discusses the origin of CF^ detection and its implications. 


2. Observations and data reduction 

We carried out an unbiased spectral survey with the IRAM 30 
m telescope on Pico Veleta, Spain, towards the source CygX- 
N63, located at (RA, Dec)j 2 ooo = (20M0"^05.2^ 41° 32' 12.0"). 
CygX-N63 is an isolated single object and can therefore be ob¬ 
served without confusion within the beam of the IRAM 30m 
(up to 30" at 3mm). The observations were obtained between 
September 2012 and January 2014 using the wobbler switching 
mode to assure good baselines, and covering so far the whole 
1, 2, and 3 mm bands (80 - 270 GHz). The EMIR receiver was 
connected to the FTS backend at 195 kHz with settings, sepa¬ 
rated by 3.89 GHz to cover the bands with a redundancy of 2. In 
order to distinguish ghosts from bright lines in the rejected side¬ 
band, we observed each frequency as pairs of settings shifted 
by 20 MHz. The data were reduced with the GILDAS packag^ 
After a zeroth-order baseline , all spectra were averaged and con¬ 
verted into main-beam temperature (Tmb) using the forward and 
main-beam efficiencies (Feff and Bgff) listed in Table [pleading 
to typical rms noise levels of 2.9, 4.2, and 8.6 mK for the 3, 2, 
and 1 mm band, respectively. 


3. Analysis and results 

3.1. Line identification and detection of CF^ (1-0) 

The analysis was performed using CASSISj^and the CDMS|^cat- 
alogue. We here report a clear detection (6.2 cr) of the J = 1^0 
line of CF^ at 102.59 GHz at vlsr = -4 km/s. The only other line 
of CF^ (J = 2^1 at 205.17 GHz) in our spectral coverage cannot 
be firmly detected because of blending with a line of CH 3 CHO 
(111 11 0 - lOi 10 0 ) (see Fig.[^. 

3.2. Total column density and excitation temperature 

Using an excitation temperature of 10 K as in [Neufeld et al.| 
( |2006| ) and |Guzman et al.| ( [2012| ), assuming local thermodynam¬ 
ical equilibrium (LTE) and optically thin lines, we derive a total 
CF^ column density of 4xl0^^cm“^ (see Appendix A for de¬ 
tails). Taking the (2^1) line tentative detection as an upper limit, 
we obtain a maximum excitation temperature of ~ 20 K. 

The line velocity of -4 km s"^ is exactly the rest velocity of 
the protostar ( [Bontemps et al.|[T010| ), which suggests an origin 
from the protostar. The study of line profiles suggests a correla¬ 
tion between the line width and the spatial origin of molecules 
(Fechtenbaum et al. in prep.). The observed full width at half 
maximum (FWHM) of 1.6 km s“^ points to an origin from the 
envelope. 

^ http:/www.iram.fr/IRAMFR/GILDAS/ 

^ http://cassis.irap.omp.eu/ 

^ http://www.astro.uni-koeln.de/cgi-bin/cdmssearch 
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Fig. 1. Spectra of the CF+ J = 1^0 and 2^1 lines on a main beam 
brightness temperature scale. The blue curves are the modeled spectra 
of the CF^ lines and the red curve is the model for CF^ and CH 3 CHO 
together. 


4. Discussion and conclusion 


4.1. First detection of CF^ towards a massive protostar 

The CF^ ion is stable and is the second fiuorine reservoir after 
HF in PDRs. Its formation route is simple (|Neufeld et al.|2005|) 
F -h H 2 ^ HF -r H and HF -r C+ ^ CF+ -r H, 
where both HF and are abundant, i.e. at the interface HI - H 2 
and more generally where ionization is strong enough to produce 
. CF^ should be an excellent proxy for in partially ionized 
regions. CF^ has so far been detected only towards two bright 
PDRs: the Orion Bar and the Horsehead nebula. The Orion Bar 
is an intense PD R partly seen edge-o n leading to a large column 
density of CF^ ( [Neufeld et al.[|2Q06 ). The Horsehead nebulae is 
a weaker PDR, but the interface is seen almost perfectly edge- 
on with a significant li mb-brightening incre asing the observed 
column density of CF^ ( [Guzman et al. 2012 1 . 

The present detection is the first dete ction towards an object 
which is most likely not a PDR (see Sect. |4.2| ). It is an interesting 
finding because it then points to the existence of a source of 
inside the CygX-N63 core. A column density of 4 x 10^^ cm“^ in 
a 25" FWHM beam corresponds to a total number of CF^ ions 
of 1.2 X 10"^^. With the abundance ratio CF^/ of the order of 
^ at hig h density and high UV field (see Fig. 4 in 


10 


et al.' 20061 , it converts to a total number of of 1 . 2 ) 


Neufeld 

_ 0^^ We 

discuss below the different possibilities to explain this amount of 
indirectly traced by CF^. 


4.2. A PDR front detected at the distance of Cygnus X? 

Since CygX-N63 is situated in the rich Cygnus X complex host¬ 
ing a large number of OB stars, the detected CF^ could poten- 
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Table 1. Observation parameters, Gaussian fits results and LTE model results. 


Line 

Erequency 

Eup/k 

^ij 

Feff 

BefF 

Beam 

Vlsr 

FWHM 

Tpeak 

Noise 


GHz 

K 




arcsec 

km 

km 

mK 

mK 

CE+(l-0) 

102.587533 

4.9235 

4.82x10-'' 

0.94 

0.79 

25.5 

-4.1 ±0.5" 

1.6 ± 0.5" 

10.5" 

2.3 

CE+(2-l) 

205.170520 

14.7703 

4.62x10-5 

0.94 

0.64 

11.4 

-4.1” 

1.6” 

<62^ 

14.7 


^ Gaussian fit results 

^ Fixed equal to the J = 1 - 0 line parameters 

^ Upper limit. LTE modeling with a 25” source size (beam size) and an excitation temperature of 10 K gives Tpeak = 24 mK. 



20‘'40”20® 40“10® 40”00" 39“50' 


a ( 2000 ) 

Fig. 2. IRAM 30m beam at 102 GHz (green circle) overlaid on a RGB 
(8, 4.5, and 3.6 jum) Spitzer IRAC image towards CygX-N63. Contours 
show the 250 jim dust emission from Herschel ( [Hennemann et al.|2012| ). 
No prominent PAH emission (red, Syum) is detected towards CygX- 
N63. 


tially originate from the PDR at the surface of the clump host¬ 
ing CygX-N63, excited by the most nearby massive stars. How¬ 
ever, CygX-N63 is located in a quiet part of Cygnus X where no 
nearby O stars can be detected (the projected distance to the cen¬ 
tral regions of the Cygnus OB2 cluster is ~20-30 pc). In compar¬ 
ison, the Orion Bar PDR is situated at less than 0.5 pc from the 
O stars of the Orion nebula cluster. To detect such a weak CF^ 
line in Cygnus X, a PDR has to be bright and should therefore 
exhibit strong PAH emission in the mid-IR. In Fig. [^showing a 
three-color Spitzer image, one sees that the CygX-N63 core ac¬ 
tually coincides with one of the darkest parts of the region with 
no sign of PAH red (8 yum) emission in the observed beam. Only 
some green (4.5 yum) weak extended features are seen precisely 
coinciding with CygX-N63. It does not trace a PDR front and is 
most probably due to outflow shock ^fro m the strong CO out¬ 
flow driven by CygX-N63 (see Sect. 4^. There is therefore no 
hint of the existence of a bright PDR at the surface of the clump 
which would explain the observed strong source of C^. 


^ There are no PAH features at 4.5 yum. Such strong excesses at 4.5 yum 
are usually referred to as EGOs (extend ed green objects) and are due to 
H2 line emission in outflow shocks (e.g. Cyganowski et al.|2008| ). 


4.3. An inner Cl I region towards CygX-N63? 


CygX-N63 is an OB star in formation. It is or will soon become 
a strong ionizing source and a strong source of C^. Interestingly 
enough, however, CygX-N63 is not yet an UC HII region. It has 
not been detected in free-free emission with the VLA. CygX- 
N63 only coincides with a 2.3 cr peak of 0.14 mJy at 8.4 GHz 
(see Appendix]^. Since carbon has a lower ionization energy 
(11.3 ev) than hydrogen (13.6 ev), a C^ region can exist before 
the development of an HII region. 

With a luminosity of 340 L©, CygX-N63 can host a B5 - B6 
main sequence star of photospheric temperature 15000 K with a 
carbon ionizing flux (betw een 11.3 and 13.6 ev) of 2.63x10"^^ 5 


( Diaz-Miller et al.[|l998 ). Assuming a r ^ law and having 44.3 
M© inside a FWHM of 2500 AU ( Duarte-Cabral et aL||2013| ), 


we get an H 2 density of 1.95 x 10^^ x(r/100AU) ^cm Using 


the recombination rate of the UMIST database ( [McElroy et al 
2013| ), assuming all carbon in atomic form and that ionization is 
dominated by carbon, we can estimate the thickness of the CII 
region (see Appendix for details). For this estimate we also 
neglect the effect of dust extinction on the ionizing radiation. 
Even with these strong assumptions the above carbon ionizing 
flux can ionize a thin layer of only 0.1 AU for an internal radius 
of 100 AU (Appendix [C]). Such a thin layer would only contain 
2.6 X 10^^ ionized carbons, i.e. ~ 500 times less than observed. 
To increase the efficiency of the ionization, the ionizing radiation 
needs to reach lower density layers where the recombination is 
less efficient so that the relatively weak ionizing radiation is able 
to keep a large amount of carbon ionized. The typical required 
H 2 density is 4 x 10^ cm“^ (see Appendix]^. Such a density is 
expected at large radii, of the order of 2000 AU, which cannot 
be easily reached by inner UV radiation because of the expected 
large dust extinction from the inner envelope. Since CygX-N63 
drives an outflow ( [Duarte-Cabral et aT]|2013| ), its envelope has 
bipolar cavities, and should have a disk-like structure in the cen¬ 
tral regions that reduces extinction in the equatorial directions. 
It is difficult, however, to expect a reduced enough extinction for 
ionizing radiation to reach 2000 AU, which is close to the phys¬ 
ical size of the envelope (see jBontemps et al.|2010| ). 


4.4. X-ray ionization from an accretion shock? 


CygX-N63 is actually a young protostar (Class 0) with a lumi¬ 
nosity dominated by accretion. An accretion shock creates a hot 
plasma which radiates half of the accretion power outward. This 
plasma has a typical temperature of 10^ K and mostly radiates 
in extreme UV, but also in X-rays. X-rays ionize gas more effi¬ 
ciently than does photospheric radiation (e.g. |Montmerle|20OT]), 
and are also less affected by extinction than UVs ( Ryter||199^ 
and could thus ionize carbon deeper in the envelope decreasing 
the extinction problem discussed in the previous section. A large 


Article number, page 3 of| 
































A&A proofs: manuscript no. aa24718-14 


fraction of X-ray ionization originates from locally emitted UV 
radiation (possibly down to carbon ionizing photons) as a result 
of local deposition of energy from photoelectric electrons gener¬ 
ated by X-rays ( [Krolik & Kallman|1983t|Glassgold et al.|2000| ). 
This has the interesting effect that X-rays, especially high-energy 
X-rays, can propagate deep in the envelope, being less affected 
by extinction than UVs, and deposit locally their energy as UV 
ionizing radiation. The total luminosity of CygX-N63 could rep¬ 
resent up to 7.2 X 10 ^® s-' ionizing photons for carbon, assuming 
that all the radiated energy is dedicated to C ionization, which is 
obviously the most optimistic view. Using the same assumptions 
as in Sect. |4.3[ we find that the carbon ionized layer for an in¬ 
ner cavity of 100 AU would then be 2.6 AU and would contain 
6.8 X 10^^ ionized carbons, i.e. ~ 20 times less than observed. 
This is sti ll a bit low, but the discrepancy is reduced compared to 
Sect. 4.3 The local density required to keep 1.2 x 10^^ carbons 
ionized would be 1.1 x 10 ^cm“^, i.e. ~ 20 times lower than 
the density at 100 AU in the pure spherical case. If the reach¬ 
able density could be a bit reduced, for instance at the surface of 
the outflow cavities, as in the irradiated outflow wall scenario of 
|Bruderer et al.| ( [2009] ), the agreement could perhaps be reached, 
providing that a significant fraction of the available energy can 
be transfered to carbon ionization as discussed above. The es¬ 
timate of the possible X-ray luminosity which could increase 


by orders of magnitude during accretion bursts (see Montmerle 
|2001| and references therein) and the overall efficiency of X-ray 
ionization at 7? ~ lOOAU accounting for X-ray absorption in 
such a complex geometry are difficult to assess precisely. We 
can only conclude that the ionizing power expected from an ac¬ 
cretion shock may reach the right order of magnitude to explain 
our CF^ detection. 


4 . 5 . Outflow shocks? 


The last possible source of ionization is outflow shocks whose 
energy is extracted from the inner regions of the collapsing en¬ 
velope. With a CO velocit y of 40 km/s and an outflo w force 
of 2.9 X 10“^ Mokm/syr“^ ( Duarte-Cabral et aL]|2013| , we get 
an outflow luminosity of 9.6 L© which could represent up to 
2 X 10"^^ s“^ carbon ionizing photons. Young jets are expected 
to radiate their energy through shocks. Three knots of 4.5 jum 
emission (shocked H 2 ) are indeed located in projection at only 
~ 3.5" fr om CygX-N63 (Fig. [^. These shocks should radi¬ 
ate in UV ( |van Kempen et al.||20()9| ). For an average projection 
angle of 57°, the shocks would be then located at a radius of 
9000 AU from the protostar where the expe cted nu density in 
the envelope is 4.8 x 10^ cm“^ (see Sect. 4^). At this density an 
ionizing flux of only 1.6 x 10^4 s-' is required to keep ionized 
1.2 X 10^^ carbons. This is less than the available ionizing power 
of 2 X 10"^^ s“^. An attenuation of the ionizing UVs in the out¬ 
flow cavities is expected, but it is reduced to only 0.24 mag for 
the expected densit}0of 3.0 x 10^ cm“^ in the cavities and for a 
distanc^of 4500 AU. The shocks in the outflow may thus pro¬ 
vide the required amount of ionizing photons to explain our CF^ 
detection. 


^ For pressure equilibrium ( [Bruderer et al.|200^ and infall and outflow 
velocities of 1 and 40 km/s respectively, the cavity density should be 
(1/40)^ of 4.8 X 10^ cm"^, i.e. 3.0 x 10^ cm"^ 

^ For an outflow opening angle of 60°, the distance to the outflow walls 
is typically of 4500 AU. For wh = 3 x 10^ cm"^, the visual extinction Ay 
is then of the order of 0.1 mag and the UV extinction of 0.24 mag using 

( [2009] ). 


the Auy/Ay ratio of 2.4 used in Bruderer et al 
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In the two possible scenarios discussed above (Sects. |4.4| and 
43] ) the produced UV radiation should ionize first hydrogen and 
create an HII region which has, however, not been detected in 
free-free emission (Fig. ED- For densities and sizes involved in 
the accretion shock case the free-free emission has to be opti¬ 
cally thick. We then estimate that the 5 cr upper limit in free-free 
of 0.3 mJy would correspond to an emitting region of ~ 20 AU 
radius (see Appendix [P]). The C ionized region is expected to be 
more extended than the HII region and given the very uncertain 
geometry and density distribution, 20 AU could be considered to 
be roughly compatible with the 100 AU scale discussed above. 
In the outflow shock scenario, for a dissociative outflow shock 
with a jet velocity of 200 km/s and = 3 x 10“^ cm“^, the free- 
free emission is optically thin with Tff = 5 x 10“"^. For an ionizing 


flux of 2 x 10 ' 


44, 


using a recombination rate of 3 x 10 "'’Hq-S 


( |Spitzer|l978l ) and for full ionization = 3x 10“^ cm“^), one 
gets EMv = 6.6 x 10^^ cm“^ leading to 0.94 mJy at 8.4 GHz (see 
Appendix[D]). This is only three times larger than the 0.3 mJy up¬ 
per limit with the VLA which could be seen as a still acceptable 
discrepancy given the large uncertainties in the above estimates. 
We can therefore conclude that in both the accretion and outflow 
shocks scenarios a CII region can be detected in CF^, while no 
free-free emission has been detected in the accompanying HII 
region. 


4 . 7 . A new probe of the early evolution of massive protostars 

In conclusion, the first detection of CF^ towards a massive proto¬ 
star may provide us with a new tool for investigating the ionizing 
power of accretion and outflow shocks. In contrast to in the 
far-IR, CF^ in the millimeter range can be observed at high spa¬ 
tial resolution with interferometer. It could offer the opportunity 
to image the development of regions even in the case of non¬ 
detection of any HII region. The indication of a possible strong 
X-ray and/or UV ionization also has important implications for 
the early chemical evolution in collapsing envelopes in the con¬ 
text of the formation of a hot core. 
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Appendix A: CASSIS LTE modeling 

CASSIS provides a Jython script which minimizes the;^^ to find 
the best parameters to adjust the observations. With a spectrum i 
of N points, the;^^ is derived as 


(Iobs,j Imodel,j) 


x- = y 

j_i rms- + cal- (Iobs,j “ Icont,j )2 


where Cbs and Imodei are respectively the observed and the mod¬ 
eled intensities, Icont is the continuum intensity, cal/ the calibra¬ 
tion uncertainty for the spectrum i and N is the total number of 
points of the spectrum. The reduced is the calculated as 


2 

Tred 


N, 


spectra 


V 

^ -^-dof’ 


Nind 


with Nspectra the number of spectra, Nind the number of indepen¬ 
dent points, and dof the degree of freedom. 

The adjustable parameters are the column density, tempera¬ 
ture, FWHM, size of the source, and vlsr- In the case of CF^, 
the column density varies between 10^^ and 10^"^ cm“^ with 500 
steps; the other parameters were fixed. The temperature was set 
to 10 K, as in |Neufeld et al.| ( |2006| ) and |Guzman et al.| ( [2012| ). The 
FWHM was measured on the spectrum as 1.6 km/s. No source 
size was assumed, so we kept the 30 m beam size of 25”. The 
velocity of the source is around -4 km/s. CASSIS also takes into 
account the beam size variation as a function of the frequency. 
The obtained is 3.22 and the column density is (4 + 1) x 
10^^ cm“^. 

The Vlsr nnd FWHM derived for the J = 1 - 0 line are then 
used to predict the intensity of the J = 2 - 1 line. We detected 
202 CH 3 CHO lines in the whole coverage and used a population 
diagram to derive a mean excitation temperature of 40 K and 
a column density of 3 x 10^^ cm“^. These values were used to 
model the CH 3 CHO lines as shown in the bottom panel of Fig. 

m 


Appendix B: VLA observations 

CygX-N63 was observed on 27 April 2003 with the 27 anten¬ 
nas of the VLA in D configuration at 8.4 GHz (band X) (project 
AB1073). The image shown in Fig. |B . 1 1 corresponds to a total 
integration time of 25 mins spread over 4.5 hr (track sharing 
technique to improve UV coverage). It has been cleaned with 
the clean algorithm of AIPS. The resulting rms is equal to 60.2 
yuJy/beam. 

In the direction of CygX-N63, a weak peak of 0.14mJy, i.e. 
at a level of 2.3 cr, is found. Below 3 cr it cannot be considered 
as a possible detection, especially because the side lobes (large 
stripes in the image of Fig. are roughly of the same intensity. 



a (2000) 

Fig. B.l. VLA image at 8.4 GHz of the same region than in Fig. with 
color scale from -0.1 to 0.7 mJy/beam. The dashed ellipse indicates the 
VLA primary beam and the green circle shows the IRAM 30m beam 
at 102 GHz. The rms in this image stays large, of the order of 0.06 
mJy/beam, due to strong sidelobes originating from to DR22, a strong 
radio source (a compact HII region) outside the imaged field (~5 arcmin 
south). 


Appendix C: Carbon ionization equiiibrium 


At equilibrium in spherical geometry from an initial radius g to 
the radius of the CII region (where mostly all carbons are ion¬ 
ized) rcii the equilibrium of carbon ionization can be expressed 
as 


2cii 




krec^c+^e-47rrdr. 


where Qcii is the total rate of ionizing photons (in s ^), ^rec 
is the carbon recombination rate that we take as equal to 
2.36 X 1 0-^^(r/300)-Q-^^exp(- 17.7/r) cm^ s'^ (UMIST 2012 
database [McElroy et al.||MT^ , and nc+ and Hq- are the densi¬ 
ties of and of electrons. For a large range of gas temper¬ 
ature (10 to 4000 K), k^Qc stays within a factor of less than 2 
around the adopted value of 1.74 x 10“^^ s“^ (ranging from 1.08 
to 2.80 X 10“^^ s“^ with the maximum obtained at ~ 60 K). 

If carbon ionization dominates and assuming C^/H equal to 
the standard C/H ratio of 1.6 x 10“"^ (if most of the carbon is in 
atomic form), we have n^- = nc+ = 1.6 x 10“"^ x nu- Assuming 
a r~^ density law and having 44 .3 Mq inside a FWHM of 2500 
AU ( [Duarte-Cabral et J^|2013| ), we here get an H 2 density of 
1.95 X 10^^ X (r/100AU)"^cm“^. For this simple model, we get 
the following expression of the extent of the ionized region in 
the inner envelope. 


^cii = 


n 

I-Tjri 


with 


Qcii 

2.85 X 1050 s-i’ 


and the radii expressed in AU. For the case of Sect. |4.3| (Qcii = 
2.63 X 10^^ s“^), one gets a very small value of 9.4 x 10“^ for rj. 
For G = 100 AU, rcn is then equal to 100.1 AU. Only the very 
first layers are ionized for values of Qcii significantly lower than 
~ 10^^ s-^ 
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The ionization equilibrium can also be expressed as 


2cii 


-f 

Jri 


^rec^e— dA^ 


with A^c+ being the number of ionized carbons. It shows that 
^rec^e- IS the fraction of ionized carbons that recombines per sec¬ 
ond. For Uq- = 1.6 X 10“"^ X we can then derive the typical 
density (assuming constant density) required to keep an amount 
of ionized carbon Nc+ ionized under an ionizing rate 2cii- Here 
is the expression of this density for our adopted numerical 
values: 

3.6x10' 

Nc+ 




.15 2cil -3 

‘ cm . 


For the case of Sect. 4.3 with 2cii = 2.63 x 10"^^ s ^ and Nc+ = 
1.2 X 10^^, one gets = 7.9 x 10^ cm“^. 


Appendix D: Centimetric free-free emission 


As discussed in Sect. |4.6[ any ionized gas is expected to emit 
free-free emission. 

The centimetric flux of free-free emission in the optically 
thin case can be expressed as a function of EMy, the emission 
measure (see [Andrej 1987 1 [Muriel et al.|1987| ), 



/ EMy \ 

i U'l 

f ^ 1 

f V 

ImJy/ 

\ 10^^ cm“^ / 

VIO^K/ ' 

18.4 GHz/ ' 

11.4 kpc / 


with 


EMy 


=///- 


dV. 


The optical depth is expressed as a function of the linear emis¬ 
sion measure EMl (= d/), 


/ £Ml \ 


f '' r * 

\ 10^^ cm“^ / 

llO^K/ ' 

18.4 GHz/ 


1.22 


which is equal to ~ 8.8 x 10 ^ for Uq- = = 2.2 x 10^cm“^, 

I = 100 AU (see Sect. |^), 10"^ K (the opacity is even larger 
for lower temperatures) and at 8.4 GHz (frequency of the VLA 
observation from Appendix]^. 

In the optically thick case the emission depends on the area 
of emission expressed in physical surface Agm if distance scaling 
is included, 



_ o fin 

( \ 

( y \2( d V 

ImJy/ 

■ ^ (100au)2 


\8.4GHz/ \ 1.4kpc/ 


)• 


which is equal to 0.3 mJy for Agm = with = 19.2 AU 
(Se ct.[4Al). _ _ 

[Curiel et~ar] ( |1987[ |1989| ) also expressed the expected flux 
and optical depth in the case of a dissociative shock as a function 
of pre-shock density and shock velocity (for > 60km/s), 

which is found equ al to 5 x 10“"^ for no = 3 x 10“^ cm“^, Vs = 
200km/s (see Sect. 4.5), 10"^ K and at 8.4 GHz. 


- 2.1 
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